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The quantum reprojection method within the standard adiabatic Born-Oppenheimer approach is 
derived for multielectron collision systems. The method takes nonvanishing asymptotic nonadiabatic 
couplings into account and distinguishes asymptotic currents in molecular state and in atomic state 
channels, leading to physically consistent and reliable results. The method is demonstrated for 
the example of low-energy inelastic Li+Na collisions, for which the conventional application of the 
standard adiabatic Born-Oppenheimer approach fails and leads to paradoxes such as infinite inelastic 
cross sections. 
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The majority of theoretical treatments of inelastic col- 
lisions involving atoms, ions, molecules, clusters, sur- 
faces, etc, is performed within the standard adiabatic 
Born-Oppenheimer (BO) approach (or simply the BO ap- 
proach), which is described, e.g., in [HE]. The approach 
is based on the separation of the electronic and nuclear 
motion. At first, the electronic fixed-nuclear Hamilto- 
nian is treated and the electronic molecular states are 
determined, then the nuclear dynamics is studied using 
the expansion of the total wave functions in terms of elec- 
tronic molecular-state wave functions. The BO approach 
gives a clear physical picture of the scattering process and 
allows one to use of the well-developed quantum-chemical 
methods and computer programs. 

Although the BO approach looks straightforward, it 
encounters severe difficulties. The problem was first rec- 
ognized in Ref. [3] and became known as the "electron 
translation (ET) problem". In fact, the name of the 
"molecular-state problem" is more adequate for this case, 
because the origin of this problem is the use of molecular 
states (which is the basis of the BO approach) in scatter- 
ing process treatments. The proposed remedies for the 
ET problem are essentially based on: (i) the inclusion of 
ET factors H] or common translation factors [HI E] into 
the expansion of wave functions, or (ii) the use of state- 
specific reaction coordinates [7HH|, (specific) hyperspher- 
ical coordinates [TO], Eckart coordinates [TTj , etc. The 
methods were reviewed in Refs. [TJ [2 EH EI] ■ The reme- 
dies lead to modifications of basis functions, potentials, 
couplings, and dynamical equations and are conception- 
ally rather complicated when compared with the original 
concept. It has been stated that it is "not possible to 
extract a meaningful scattering matrix" |12j . and finally 
the problems have been interpreted as conceptional lim- 
itations of the entire BO approach [2], Nowadays, this 
fundamental problem still presents unresolved features, 
e.g., infinite scattering lengths in ultralow-energy colli- 
sions [14 . Efforts to solve the ET problem have been 
continued, and, in particular, the one-electron quantum 
reprojection method has been derived |15H17| . which is 
conceptionally rather simple and uses BO molecular po- 
tentials and couplings. This method is generalized in the 



present paper for a multielectron case. 

In most applications, however, the ET problem is sim- 
ply neglected, that is, (i) all asymptotic nonadiabatic 
couplings are cut off at a finite internuclear distance, and 
(ii) the asymptotic boundary conditions are taken in the 
BO (see below) or similar form, see, e.g., [TB]- Let us 
refer to this application of the standard adiabatic BO 
approach as the "conventional BO method", because the 
majority of cross-section calculations is carried out with 
this procedure, in contrast to the above mentioned meth- 
ods which take the ET problem into account. When ap- 
plying the conventional BO method, it is assumed that 
the approximations described above give negligible er- 
rors, at least for low collision energies, see, e.g., [T51 fTS] . 
Nevertheless, as shown in Refs. [I7J[20], even at low col- 
lision energies the conventional BO method leads to the 
paradoxes such as nonzero nonadiabatic transition prob- 
abilities and nonzero inelastic cross sections for noninter- 
acting collisional systems, e.g., n + H. 

Taking as an example Li+Na collisions, the present 
paper shows that in some cases no proper cutoff can be 
found to get reliable results and that the above mentioned 
approximations can lead to errors which are several or- 
ders of magnitude larger than correct values even at low 
collision energies. Moreover, it is shown that in case of 
nonzero asymptotic couplings, which are the rule rather 
than the exception of the BO approach, the conventional 
BO method leads to infinitely large inelastic cross sec- 
tions. Therefore, the coupling-cutoff procedure can give 
any value for an inelastic cross section. 

For the sake of simplicity, let us treat atomic collisions 
in E molecular states. Within the BO approach the total 
wave function ^ j m., (r, R) is expanded as 

9j Mj (r, R) = Yj Mj (Q, $) <f>k (r, R) , (I) 

k 

</>fc(r, R) being the electronic molecular-state wave func- 
tions, r and R being the sets of electronic and nuclear 
coordinates, J, Mj being the total angular momentum 
quantum numbers. This results in a system of coupled 
channel equations (CCE) for radial nuclear wave func- 
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tions Fk(R), see, e.g., pQ. Nonadiabatic transition prob- 
abilities are then calculated in the asymptotic (i? — > oo) 
region, where the conventional BO method assumes that 
an incoming/outgoing current in a single atomic state 
proceeds completely into a single molecular state and vice 
versa. The asymptotic boundary conditions for the total 
wave function read 
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with the wave numbers Kj and the outgoing/incoming 
amplitudes in the atomic-state channel j. The con- 
ventional BO method assumes the following outgoing and 
incoming asymptotic (R — > oo) BO wave functions 



(3) 



The CCE are solved numerically from zero up to an up- 
per integration limit Rq, which is large enough to cal- 
culate transition probabilities based on the asymptotic 
wave functions (|3| PQ. 

The use of electronic molecular states leads to the fol- 
lowing fundamental features. Both radial and rotational 
nonadiabatic couplings, see, e.g., [H [15j EU E2] , as well 
as the form of CCE jTS], [55] depend on the origin of the 
electron coordinates, although the CCE themselves are 
independent on the origin choice Q]3 [22123]. The CCE 
take their standard and simplest form in Jacobi coordi- 
nates, where the electrons are measured from the center 
of nuclear mass (CNM) (neglecting the mass-polarization 
term) . The asymptotic values of the radial nonadiabatic 
couplings calculated with the electron origin at the CNM 
read, see, e.g., [TSllTo] . 
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01^|fc)oo = Ik « M(oo) - V k (oo)] (j\df\k) , (4) 



(j\dg \k) being the atomic transition dipole moment, m 
being the electron- nuclei reduced mass, Vj(R) being an 
adiabatic potential. The scalar factors 7^ depend on 
which nucleus an active electron is bound in the asymp- 
totic region: 
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an electron bound with A 
an electron bound with B. 



(5) 

A and B label of nuclei with masses Ma and M b ■ About 
rotational couplings, see |16j . It is seen that some radial 
nonadiabatic couplings remain nonzero in the asymptotic 
region even for the noninteracting model system n+H |17l 
[20] 122] . It is worth to emphasize that choosing another 
electron origin (e.g., at one of the nuclei) does not help to 
avoid nonzero asymptotic couplings in the CCE, because 
the same nonzero values appear in the equations due to 
new terms in the Hamiltonian [TSJ [22] . According to 
the BO approach, nonzero couplings provide transitions 
between molecular states even at R — > 00. 



In fact, the nonvanishing asymptotic couplings are a 
consenquence of a more fundamental shortcoming. The 
coordinates used to describe molecular states of the colli- 
sion complex at small and intermediate distances are not 
suited for the description of the free atoms in the asymp- 
totic region. The correct asymptotic incoming/outgoing 
wave functions Q] 1 E3 HU [H] 



ex.p(±iK j R° 
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are written in another set of Jacobi coordinates and dif- 
ferent from the BO functions (J3J. The vector R|* con- 
nects the centers of mass of the atoms, in contrast to R 
which connects the nuclei. 

The reprojection method for multielectron collision 
systems consists in the following. The vector R"' can 
be written as follows 



Rf = R + bj , 



(7) 



where the vector hj is equal to a sum over all electrons 
'>, w (r Q -7?R) , (8) 



index a labeling the electrons, r = {r Q }. m" is the 
electron- nucleus reduced mass in the channel j: m" = 
m e MA/(m e + Ma), if the electron a is bound to nu- 
cleus A, and m" = m e M B /(m e + M B ), if it is bound 
to nucleus B; m e being the electron mass; M being the 
reduced mass of the nuclei. The shift bj depends on the 
asymptotic electron rearrangement in the channel j and, 
hence, can be different for different channels. It is small, 
but it does not vanish at infinity and therefore should be 
taken into account. 

A single term in the expansion (JlJ , describing a current 
in the molecular state k, does not coincide with a single 
term in Eq. ([2]), describing a current in the correspond- 
ing atomic state, see Eq. Q. An incoming/outgoing cur- 
rent in a single atomic state is distributed among several 
molecular states and vice versa. Projecting of the atomic- 
channel asymptotic wave functions ([6| on the molecular 
asymptotic wave functions ^ gives 



± = cxp{±iK j R) 
1 R 
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where the elements of the matrices ty^ represent the re- 
projection coefficients 



hj = ( k \ ex P( ±iK j b j z )\j) c 



(10) 



At low collision energies these matrix elements can be 
approximately evaluated via corresponding atomic tran- 
sition dipole moments and furthermore, taking into ac- 
count Eq. Q, via asymptotic values of the derivative 
couplings calculated in the Jacobi molecular coordinates 
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FIG. 1: The adiabatic potentials (upper panel, a) for the 
three lowest LiNa( 1 E + ) states (X, A, C) and the radial nona- 
diabatic couplings (low panel, b) between these states. 
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FIG. 2: The transition probabilities for Li(2s -> 2p) +Na(3s) 
(upper panels) and Li(2s) + Na(3s — > 3p) (low panels) exci- 
tation calculated for E = 5 eV and J = by means of the 
conventional BO method (the solid black and red lines) and 
by means of the reprojection method (the dashed blue and 
green lines) as a function of the upper integration limit Rq. 
The right panel shows the probabilities in the enlarged scale. 



where all values are taken in the asymptotic region. So, 
the asymptotic couplings are responsible for the correct 
asymptotic wave functions in the r, R coordinates. 

A numerical solution of the CCE for a molecular chan- 
nel at R — ¥ oo gives a superposition of the BO asymptotic 
wave functions in atomic channels. Within the reprojec- 
tion method, the CCE with nonzero asymptotic nona- 
diabatic couplings are integrated numerically from zero 
up to a large distance Rq resulting in the i?-matrix (R). 
Finally, taking into account Eq. Q, the S'-matrix is ex- 
pressed via the i?-matrix as follows 

S = (-l) J exp(-iiLRo)^r 1/2 (f+iRfK) 

x (|+ - igt+g)^ 1 K 1/2 exp(-iJCRo) ■ (12) 



The formula ( 12 ) is valid for multielectron collisions and 



turns into the formula derived previously in the limiting 
one-electron case [T51 ITB] . The difference between the 
conventional BO method and the reprojection method is 
in the presence of the ^-matrices in the latter instead 
of the unit matrix in the former, see Eq. (12 1. The im- 



plementation of the reprojection method (also called as 
the t-matrix method) is not more complicated than the 
implementation of the conventional BO method and does 
not require additional input data if Eq. (11) is used. 



Let us consider Li+Na collisions. The adiabatic poten- 
tials for the three lowest LiNa( 1 S + ) states (X, A, and C) 
and the radial nonadiabatic couplings between them are 
plotted in Fig. [T] The ab initio potentials are taken from 
Ref. [21] , while the nonadiabatic couplings have been cal- 
culated in the present work by means of the MOLPRO 
package (the ab initio MRCI method). Fig. [I] clearly 
shows that two nonadiabatic couplings remain nonzero in 
the asymptotic region in agreement with Eq. Q. Their 



values are not negligible compared with the typical maxi- 
mum value of w 0.2 a.u. Moreover, the X-A coupling has 
its maximum value in the asymptotic region (except for 
the range R < 2 a.u., which is not important for tran- 
sitions). Thus, no proper cutoff can be found to get a 
reliable result. 

The transition probabilities Pif(J,E) for the Li(2s) + 
Na(3s) -5> Li(2p)+Na(3s) and Li(2s)+Na(3s) ->■ Li(2s) + 
Na(3p) excitation processes are shown in Fig. [2] for the 
collision energy E — 5 eV and J = as a function of 
the upper intergation limit i?o- It is seen that Pif(J,E) 
calculated by means of the conventional BO method os- 
cillate between roughly zero and relatively large values 
with increasing Rq. The variations represent nonadia- 
batic transitions between molecular states at large R and 
are the consequence of the nonzero asymptotic couplings. 

The reprojection t-matrix method yields transition 
probabilities which are independent on the upper inte- 
gration limit, when it is large enough, sec Fig. [2] The 
transition probabilities are several orders of magnitude 
smaller those obtained by the conventional BO method 
and are not equal to averaged values of the latter. Nona- 
diabatic transitions between molecular states still remain 
at an arbitrary large R, but they do not produce tran- 
sitions between atomic states at large distances. The 
t-matrices correct the S'-matrix (12 1. Thus, transitions 



between atomic states in the asymptotic region are un- 
physical, while transitions between molecular states in 
the same region are physical and represent a fundamen- 
tal feature of the BO approach. 

The same transition probabilities Pif(J,E) as a func- 
tion of J are plotted in Fig.[3]for E = 5 eV, Rq = 500 a.u. 
The conventional BO method gives probabilities which 
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FIG. 3: The transition probabilities for Li(2s -> 2p) +Na(3s) 
(a) and Li(2s) + Na(3s — > 3p) (b) excitation as a function of 
the total angular momentum quantum number J for E = 
5 eV. The thin black and red lines show the probabilities 
obtained by means of the conventional BO method; the thick 
blue and green lines depict the probabilities calculated by the 
reprojection method and multiplied by 1000 (a) and 200 (b). 
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FIG. 4: The Li(2s -> 2p) + Na(3s) (black and blue lines) 
and Li(2s) + Na(3s —¥ 3p) (red and green lines) excitation 
cross sections calculated by means of the conventional BO 
method (dashed lines) and by means of the reprojection t- 
matrix method (solid lines) as a function of the maximum 
value of the angular momentum quantum number J m ax for 
E = 5 eV. 

remain oscillating with increasing J due to the nonzero 
asymptotic couplings: at any J a centrifugal term does 
not prevent reaching a nonadiabatic region. In con- 
trast to this, the transition probabilities obtained by 
means of the reprojection method are substantial only 
within a limited range of J, roughly up to J ~ 500 for 
E = 5 cV. Note the probabilities calculated by the re- 
projection method and plotted in Fig. [3] are multiplied 
by 1000 and 200. 

The inelastic cross sections are calculated as a sum 



over J from till infinity or a value J ma x where a con- 
vergence is reached. If the range of J with nonzero tran- 
sition probabilities is unlimited, as shown in Fig. [3] for 
the conventional BO method results, convergence cannot 
be reached and cross sections infinitely increase with in- 
creasing of the upper summation limit J m ax, as depicted 
in Fig. [4] Thus, nonzero asymptotic couplings lead to 
infinite inelastic cross sections obtained with the conven- 
tional BO method. This conclusion holds even in case of 
small nonzero asymptotic couplings. 

In the framework of the reprojection t-matrix method, 
the nonadiabatic transition probabilities Pif(J,E) are 
nonzero only within a limited range of J (as it must be) , 
which leads to the convergence of the cross sections with 
increasing of J m ax, see Fig.|4j and finally to the finite val- 
ues of the inelastic cross sections. All remedies for the ET 
problem are supposed to lead to the convergence. Thus, 
the conventional BO method has its limitation both in 
the formalism and in its applications, while the standard 
adiabatic BO approach with any ET remedy is free from 
such limitations. 

It should be mentioned that a lack of convergence with 
respect to an integration limit variation within the con- 
ventional BO method at high energies was noticed in 
Ref. [5] for another processes. In those papers it was 
found that some variation of integration limits results in 
increasing of the total cross sections up to a factor of 3 at 
high collision energies, while at low energies (< 300 eV) 
"the total cross sections without common translation fac- 
tors are only slightly different." In the same papers it 
was found that the inclusion of the common translation 
factors leads to the convergence with respect to an inte- 
gration limit variation. In the present work it is found 
that within the conventional BO method there is no con- 
vergence with respect to both an upper integration limit 
R variation (Fig. [2]) and a variation of the maximum 
value of the angular momentum quantum number J max 
(Fig. |4| even at low collision energies, while the repro- 
jection t-matrix method provides the convergence with 
respect to variation of both R and J m ax, see Figs.[2]and 
[4] Moreover, it is shown that at low collision energies, 
the results of calculations with and without taking into 
account the ET (or molecular-state) problem differ not 
only by a few times, but by several orders of magnitude. 

It has thus been demonstrated that the conventional 
BO method applied to collision processes with nonzero 
asymptotic nonadiabatic couplings, which are fundamen- 
tal features of the BO approach, leads to paradoxes such 
as infinite inelastic cross sections even at low collision 
energies. The reprojection method takes into account 
nonzero asymptotic couplings, which are responsible for 
correct asymptotic wave functions, and distinguishes the 
asymptotic currents in molecular and atomic state chan- 
nels, providing physically consistent and reliable results. 

The author thanks Prof. W. Domcke for fruitful dis- 
cussions. The work was supported by the Russian Foun- 
dation for Basic Research (Grant No. 10-03-00807-a) and 
the University of Rome "La Sapienza". 



[1] A. Maci'as and A. Riera, Phys. Rep. 90, 299 (1982). 

[2] B. H. Bransden and M. R. C. McDowell, Charge Ex- 
change and the Theory of Ion - Atom Collisions. (Claren- 
don, Oxford, 1992). 

[3] D. R. Bates and R. McCarroll, Proc. Roy. Soc. London, 
Ser. A 245, 175 (1958). 

[4] W.R. Thorson and J.B. Delos, Phys. Rev. 18, 117 (1978). 

[5] S. B. Schneiderman and A. Russek, Phys. Rev. 181, 311 
(1969); L. F. Errca, L. Mendez, and A. Riera, J. Phys. B 
15, 101 (1982). 

[6] R. J. Allan and J. Hanssen, J. Phys. B 18, 1981 (1985); 
R. J. Allan, A. Bahring and J. Hanssen, ibid, 18, 1999 
(1985). 

[7] M. H. Mittleman, Phys. Rev. 188, 221 (1969). 

[8] W. R. Thorson and J. B. Delos, Phys. Rev. 18, 135 

(1978); J. B. Delos, Rev. Mod. Phys. 53, 287 (1981). 
[9] M. Gargaud, R. McCarroll, and P. Valiron, J. Phys. B 

20, 1555 (1987); D. Rabli and R. McCarroll, J. Phys. 

B 38, 3311 (2005); M. Roudjane, R. McCarroll and D. 

Rabli, J. Phys. B 40, 2491 (2007); P. Barragan et al, 

Phys. Rev. A 81, 062712 (2010). 
[10] E. A. Solov'ev and S. I. Vinitsky, J. Phys. B 18, L557 

(1985); J. Macek et al, Phys. Rev. A 35, 3940 (1987). 
[11] J. Robert and J. Baudon, J. Phys. B 19, 171 (1986). 
[12] R. McCarroll and D. S. F. Crothers, Adv. At. Mol. Opt. 

Phys. 32, 253 (1994). 
[13] L. F. Errea et al, J. Phys. B 27, 3603 (1994). 
[14] L. Wolniewicz, Phys. Rev. A 68, 042717 (2003). 
[15] J. Grosser, T. Menzel, and A. K. Belyaev, Phys. Rev. A 

59, 1309 (1999). 
[16] A. K. Belyaev et al, Phys. Rev. A 64, 052701 (2001). 
[17] A. K. Belyaev, Physica Scripta 80, 048113 (2009). 
[18] M. Stenrup, A. Larson, and N. Elander, Phys. Rev. A 

79, 012713 (2009). 
[19] C. H. Liu et al, Phys. Rev. A 81, 012707 (2010); J. Loreau 

et al, Phys. Rev. A 82, 012708 (2010). 
[20] A. K. Belyaev, D. V. Vlasov, and A. M. Kas'yanova, Opt. 

Spectrosc. 103, 920 (2007). 
[21] T. G. Heil, S. E. Butler, and A. Dalgarno, Phys. Rev. A 

23, 1100 (1981); L. F. Errea et al, J. Chem. Phys. 84, 

147 (1986); B. Zygelman and A. Dalgarno, Phys. Rev. A 

33, 3853 (1986). 
[22] A. K. Belyaev, A. Dalgarno, and R. McCarroll, J. Chem. 

Phys. 116, 5395 (2002). 
[23] B. Zygelman et al, Phys. Rev. A 46, 3846 (1992); G. J. 

Bottrell, Nucl. Instr. Meth. Phys. Res. B 79, 173 (1993). 
[24] I. D. Petsalakis, D. Tzeli, and G. Theodorakopoulos, J. 

Chem. Phys. 129, 054306 (2008). 



